Displacement Profile of Charge Density Waves and Domain Walls at Critical 

Depinning 
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The influence of a strong surface potential on the critical depinning of an elastic system driven 
in a random medium is considered. If the surface potential prevents depinning completely the 
elastic system shows a parabolic displacement profile. Its curvature C exhibits at zero temperature 
a pronounced rhombic hysteresis curve of width 2/ c with the bulk depinning threshold f c . The 
hysteresis disappears at non-zero temperatures if the driving force is changed adiabatically. If the 
surface depins by the applied force or thermal creep, C is reduced with increasing velocity. The 
results apply, e.g., to driven magnetic domain walls, flux-line lattices and charge-density waves. 

PACS numbers: 75.60.-d, 74.60.Ge 



The driven viscous motion of an interface in a medium 
with random pinning forces is one of the paradigms of 
condensed matter physics P, This problem arises, 
e.g., in the domain wall motion of magnetically or struc- 
turally ordered systems with impurities 3] or when an in- 
terface between two immiscible fluids is pushed through 
a porous medium Closely related problems are the 
motion of other elastic systems like a vortex line in an 
impure superconductor 5] , of a dislocation line in a solid 
or driven charge density waves (CDWs) Q. For a con- 
stant external driving force this problem has been con- 
sidered close to the zero temperature critical depinning 
threshold [3, H UiJ liil an d in the creep region 0, 02 • 

It was a tacit assumption of these investigations that 
the motion of the elastic system is not hindered by effects 
from surfaces or internal grain boundaries. Surface bar- 
riers are however known to be relevant in all cases men- 
tioned above. In superconductors they prevent the pen- 
etration of new flux lines into the probe In CDWs 
normal electrons have to be converted into those con- 
densed in the CDW by a phase-slip mechanism which 
is essentially a nucleation process 0, . The motion 
of domain walls may be hindered by a variation of the 
width of the sample such that position of minimal width 
are preferred etc. 

It is the aim of the present Letter to consider the 
effect of a strong surface pinning potential in addition 
to the weak bulk random pinning. It turns out that, 
starting from a flat interface, at T — and increas- 
ing the driving force / to / > f c the average displace- 
ment profile shows a parabolic profile with a mean cur- 
vature C(f) = (/ — / c )/r. In more general situations 
C(/, i) exhibits a pronounced hysteretic behavior. At 
non-zero temperatures C(/, t) increases with time and 
reaches asymptotically its value of the pure system //T. 
We further determine the reduction of the curvature in 
the case that the surface is depinned due to a sufficiently 
large driving force or due to thermally activated processes 
at the surface. The latter mimic also phase slip processes 
in CDWs. 



Model and zero temperature critical depinning. — We 
focus on a simple realization of the problem. The equa- 
tion of motion of a D-dimcnsional field y(x, t) describing 
an interface profile in the case of domain walls or a phase 
profile in the case of CDWs is given by = 
where TL denotes the Hamiltonian of the system: 
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7 and T denote the mobility and the stiffness constant 
of the elastic object, respectively, and / is the driving 
force which is assumed to change only adiabatically. The 
potential includes a random force and a surface contri- 
bution 



d<p'g( X , P (x)] + V a (^)p(x) . (2) 



The random force — <jf(x, tp) is assumed to be Gaussian 
distributed with (g) d — and (<?(x, <p)g(x', <p')) d — 
<j(' D )(x — x')Ao(<y3 — if') where (. . . } d denotes the ran- 
dom average. For domain walls Ao(y>) = Ao(— ip) is an 
analytical monotonically decreasing function of ip which 
decays to zero over a finite distance /. For CDWs 
g oc sm(ip — a(x)) with a random phase a(x) S [0, 2n[ and 
therefore Ao(<p) is periodic with Aq(<^) = Ao(ip + 27rZ). 

The surface potential V s ((p) is assumed to act only 
in the vicinity (a <C L) of x\ = and x± — L, e.g., 
p(x) = e~ Xl l a + e( Xl ~ L ^ a , and favors the values of 
ip(0,xj_) and <p(L,x±) at 27rZ. The details of the inter- 
action between the elastic system and the surface depend 
on the specific system under consideration. We will re- 
strict ourselves here to a periodic surface potential which 
has applications in type-II superconductors and may also 
serve as a first step for the treatment of conversion phe- 
nomena in CDWs. 

The case V s = was considered previously in 0, |j| 0, 
Hlj |. It was shown that at zero temperature the system 
undergoes a depinning transition at a critical value f c . 
For f > f c the velocity v = (ip) increases as v ~ (f — fc) 13 



2 



with the critical exponent (3 calculated in an expansion in 
D = 4 — e dimensions. The average displacement profile 
is macroscopically flat. At non-zero temperatures the 
dcpinning transition is smeared out and goes over into a 
creep motion for / -C f c 0- 

In this Letter we will consider the opposite case where 
a strong surface potential V s slows down or prevents 
completely the motion of the elastic object. Then the 
displacement profile becomes parabolic with a history- 
dependent curvature C(f,t) . The steady state solution 
for the average phase is given by 



( l p)=vt + ^-{L-x 1 )x 1 



(3) 



where C s (f) = C s (f, t — > oo) is the saturation value of 
the curvature. 

Infinite surface barriers. — We begin with the case 
V s — > oo, where the depinning transition is depressed. 
To determine C{f,f) we first apply perturbation theory 
Using the decomposition p(x) = po{x) + < Pi( x ) with 
(ipi(x)) d = in the equation of motion and expanding 
<7(x, ipo + ipi) to linear order in ipi we get after averaging 
over the disorder 

i^o = -rC(t)+f+{g ip {x,<p {x,t))cp 1 (x,t)) d , (4) 

where g v (x,p) — J=^g(x,p). Calculating ipi also to first 
order of g we get from © and (@J C s = Co = f /T since 
Aq(0) = 0, i.e., there seems to be no influence of the dis- 
order. Here, the situation is completely analogous to that 
at the conventional depinning transition 8] . However, as 
we know from critical depinning, this is the situation be- 
low the Larkin scale L p . 

Next we discuss renormalized perturbation theory 
starting from a situation where C s = 0. As long as / < f c , 
the elastic object is pinned and boundary pinning does 
not matter, hence C(f, t) = 0. At / = f c the elastic 
object is in the same critical state as at the depinning 
transition. Therefore we can use the results of the pre- 
vious renormalization group calculation in this case. As 
a result 7 and Aq(z) are replaced there by the renormal- 
ized, momentum dependent quantities 

7 (p) ~ 7(pLp)- 2+z , (5a) 
A» « K D \TIIL^- D -^^ (<p(pL P )</l)) (5b) 

C denotes the roughness exponent which was calculated 
for domain walls to order e = 4 — D in 0, 13 and recently 
to 0(e 2 ) For CDWs ( = 00). The most important 
feature of A p (ip) is that A* (ip) has a cusp-like singularity 
at the origin. The renormalized equation for C(f,t — » 00) 
is given by 

00 1 

rc(/, t) = f + f p uj p J d£ Jdpx (6) 


x -i +z -c e -^r*' A *' (j_ c(t) + c{t _ t ')]^l( Xl - l)) , 



where p = pL p , tu p = jf p /l and f p = ITL~ 2 . After 
having increased / adiabatically to a fixed value slightly 
larger than / c , C(/, i) saturates for t — > 00 and hence 
the difference C(t) — C(t — t 1 ) vanishes. As a result the 
argument of A*' also vanishes and the right hand side of 
© becomes independent of x\. Since C(t) > C(t — t') the 
argument of A*' approaches zero from positive values. 
Thus we get for the saturation value C 8 (f) 



Cs(f) = 



f-fc I f-fc 



L v f P 



fc = 



fp 



2-C 



A*'(0+^ 



(7) 

One can understand this result in the following way: Us- 
ing the decomposition ip = tpo + ipi in the asymptotic 
region, where C(t) saturates, the equation of motion can 
be written as 



1 



^ =rvVi + /-rc + Sl (x, w ) 



(8) 



where gi(x,(fi) = g(x, ipo(x) + Pi(x)). gi(x,cp) and 
g(x, ip) have the same statistical properties. According 
to © the force acting on the field p\ is now reduced by 
the curvature force — TC. The depinning of the </?i-field 
seems hence to occur at / f f c = f c + TC S . However, 
since the boundary conditions fix ipi(0) — <Pi(L) = and 
hence {ipi) = for all values of /, the system is always 
at its depinning transition, which implies Q. Starting 
from some f < f c and C s = 0, C s will stay at this value 
until / reaches f c . For / > f C) C s obeys 0. The same 
argument can be used for negative forces / < 0. Then 
we find for / < -f c : TC s {f) =f + f c = -(\f\- f c ) since 
A*'(0-) = -A*'(0+). 

A scaling argument supports the validity of eq. Q 
to all orders in g: Close to the V s = depinning 
transition the correlation length £ diverges as £ ~ 
Lp ((/ — fc)/ fc) V ■ For L' < £ the roughness - the mean 
square displacement of a piece of linear size V of the elas- 
tic object - scales as w 2 (L') « l 2 (L'/L p ) 2( HUES EH. 
If we choose the system size L m £ we expect that the 
roughness scales as the height of the parabolic ip— profile 
on the same scale, u>(£) « C£ 2 , which is indeed mlfilled 
if we use the scaling law v = 1/(2 — C) 
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(9) 



Hysteresis. — Next we consider the case that we in- 
crease / adiabatically from / < f c to a value / max , where 
C(f,t) reaches C max , and then decrease f again. In this 
case C(/, t) < C(f,t — t') and hence the argument of A*' 
becomes negative. Instead of Q we get from (JSJ 



FCmax — /max fc / ~t~ fc 



(10) 



The effective force acting on the elastic object is now 
given by / — rC max . Further decreasing /, there is 
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no change of C(f,t) until the effective force reaches the 
threshold - f c = f- TC max = f - (/ max _- / c ). According 
to the last relation this happens at / = / max = f m ax~ %fc- 
Analogous arguments can be used for reversing the fields 
from / < to / > 0. Thus C s undergoes a hystere- 
sis which consists of the two parallel segments given by 
C S T = (/ T fc) and two horizontal segments determined 
by C max = (/ max - f c )/T and C min = (/min + /c)/I\ re- 
spectively. Indeed, similar hysteresis effects of the strain 
have been observed in CDWs [lfij . 

These findings are fully supported by numerical sim- 
ulations, as shown in fig. ^ For integrating the equa- 
tion of motion, the x-coordinate is discretized with a 
lattice constant a and the simulation time is measured 
in units of a time To (the dimcnsionless lattice laplacian 
for D = 1 is given by 'V 2 fi = fi+i + <fi—i — 2<pi, with 
lattice sites i = 0, ...,L). a and tq are chosen such 
that = 1 and the dimensionless stochastic forces 

To7p(x, ip) E [—1/2,1/2] (the dimensionless driving force 
is T07/). 




f 



FIG. 1: Hysteresis of C(f) at T = for an one- dimensional 
interface. The driving force is first increased to / max = 0.6 
or /max = 0.8, respectively, and then decreased to —f c £s 
—0.27. The arrows show the direction of the hysteresis. The 
numerical simulation was done for an interface with length 
L — 1000 and averaged over 300 disorder configurations. 

Curvature at finite temperature. — Next we want to 
consider the problem of finite temperatures. Changing 
/ only adiabatically we may use equilibrium statisti- 
cal mechanics. It is convenient to go over to the field 
</?(x) = (ys(x) + JpXi(xi — L). The Hamiltonian rewrit- 
ten in if has the same statistical properties as the initial 
one (QJ, since Vr(x, if) — — J v d<p'g(x, if') is a random 
function of both arguments. This can most easily seen 
by using the replica method 0. The disorder averaged 
free enthalpy follows from the replica Hamiltonian 

H " = \ E / \(V£a) 2 S a , b - ^R&a - <h)\ , (11) 
a,b=l ^ x ^ ' J 



with {V R ( X ,f)V R (x' : if')) d = <5< D )(x - x')R{<p - iff) Ap- 
parently, the replica Hamiltonian is the same as that fol- 
lowing from Q. It is worth to mention that this is true 
only if the random potential Vr(x, if) is strictly uncorre- 
lated in x-direction. The application of surface barriers 

1 /2 

implies therefore C = f/T and (((^(x) - ^(x')) 2 )^) ' ~ 

I (L/Lp) where £ denotes the equilibrium roughness ex- 
ponent corresponding to Hamiltonian I jllj l . Thus the dis- 
placement profile is the same as in the pure case. For 
non-adiabatic changes of /, traces of the T = hystere- 
sis are expected to be seen at non-zero temperatures (cf. 

fig. El- 
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FIG. 2: Simulation-time resolved coefficient C(t) for a driving 
force of / = 0.2 at various temperatures. The simulation was 
done for a system of length L — 1000 and for one disorder 
configuration of CDW-type for each temperature (see text). 

The numerical solution of equation of motion with 
thermal noise at finite temperatures and V s = oo is in 
agreement with these analytical considerations. Fig. [21 
shows the coefficient C(t) as it approaches its saturation 
value C s = f I (2r) with time. Strictly speaking, we are 
not in a steady state until C(t) has reached its saturation 
value and hence the phase profile deviates slightly from 
the parabolic shape. In fig. [3 C(t) is the least square fit 
to the profile. Note, that for low temperatures (T < 5.0 
in the simulation, where T is the dimensionless variance 
of the thermal noise) this approach is very slow, notice- 
able by the occurring steps, triggered by avalanches, even 
at large times. For high T (T = 5.0), one sees that C(t) 
fluctuates around the saturation value due to thermal 
noise. Therefore the T = hysteresis of C vanishes at 
finite temperatures. 

Critical Depinning. — So far the surface potential was 
assumed to fix the value of the displacement field ip at the 
surfaces x\ = and x\ — L. We will now assume that 
the surface potential is reduced (or temperature is raised 
from zero) such that a macroscopic motion of the elastic 
object is possible. To determine the mutual interaction 
between the bulk and the surface we have to consider 
the effective equation of motion of the surface. Denoting 
if(0,x±) — <y9 s (xj_) the effective equation of motion of 
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the surface field can be written as (a is assumed to be of 
the order of the lattice spacing) 



i 



-<Pa 



TC- 



2a 



(12) 



An analogous equation can be written for ip(L, xj_). In 
(|12|1 we have replaced the force resulting from the dis- 
placement in the bulk by the corresponding average force. 
In the steady state V^9? s = and eq. ljT2|) has a depin- 
ning threshold f SiC ^> f c determined by 



fs.c jmaxV/^) = 



(13) 



For / > f sc > f c the macroscopic velocity is given by 
the steady state solution v — ip s which follows from inte- 
grating l|12fl with V^(y3 s = 0. The corresponding solution 



v(t) 



v p $ 



YC^f + fa 2 
TV 



(14) 



depends of course on the specific form of the surface po- 
tential, = 7 f p . Ea. l|14l) has to be combined with the 
effective equation for the bulk (/ > f c ) || 



Vn 



1//3 f _ f L 2 

jp » 



(15) 



which follows from (|5"a|) and (J7J). Eqs. (jl*4"|) and l(TB|) de- 
termine both the velocity and the curvature C as a func- 
tion of the driving force. If we increase / from / = 
with C — 0, C remains zero until we reach f c . For 
fc < f < fs,c, C obeys Q. At f StC the elastic object 
is depinned and curvature is reduced with increasing ve- 
locity, according to (|T5|) . If the surface potential is pe- 
riodic, also v(t) will be periodic and the bulk depinning 
transition is slightly smeared out 01- We will assume 
that this effect is weak. In principal it can be avoided by 
adding some randomness to the surface potential. 

Nucleation and Creep. — At finite but low tempera- 
tures the surface field may exhibit a creep motion even if 
/ «C f s ,c- Creep proceeds via the formation of droplets at 
the surfaces x\ = and x\ = L, inside which cp is changed 
by 2-7T, respectively, with respect to the bulk value of tp. 
The droplet consist of a cylindrical piece (the cylinder 
axis is perpendicular to x\ = 0) in the surface layer of 
height a and radius R and an attached semi-sphere with 
the same radius. The width of the droplet wall confining 
the cylinder is of the order a' — a/^/V". 

Keeping only the leading order 
get for the energy of the droplet 
2TR°- 2 { ^/Vj 7 + In § - nCRL - fR 2 /T} 
droplet size Re <C L follows as R c w 
R c = (CL)^ 1 in D = 3 or D = 2, respectively. In 
deriving Ed P , we have neglected the contribution from 



terms we 
Ed P {R) = 
The critical 

tW/tCL) or 



the disorder which is correct as long as R c < L p , i.e., 
L > L p . The nucleation rate of droplets and hence the 
creep velocity is given by 



VD=2 



= A' 



CLa' 



bT/t 



(16) 

which replaces ifTHl in the case / <C f s ,c, T > (cf. [lj]). 
The present treatment is too crude to give the coefficients 
A,A',B,B'. Again, l|ltj|) has to be considered together 
with l|15fl to determine C and v. In CDWs, where tp can 
be multi-valued, nucleation processes also occur deep in 
the bulk ^ij. The droplet energy then does not con- 
tain a term ~ /, leaving the relations Q16f) essentially 
unchanged. 

To conclude we have shown that surface pinning of 
impure elastic systems lead to an onset of curvature C 
only above a threshold value f c of the random force. In 
general, C exhibits a pronounced hysteresis. The curva- 
ture is reduced above the surface depinning transition or 
at finite temperatures when nucleation processes at the 
surface allow for creep motion. 
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